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ON THE RIEMANN-HILBERT FACTORIZATION PROBLEM FOR 
POSITIVE DEFINITE FUNCTIONS 

DAN KUCEROVSKY, AMIR T. P. NAJAFABADI, AND AYDIN SARRAF 


Abstract. We give several general theorems concerning positive definite solutions of Riemann- 
Hilbert problems on the real line. Furthermore, as an example, we apply our theory to the 
characteristic function of a class of Levy processes and we find the distribution of their 
extrema at a given stopping time. 


1. Introduction 

The Wiener-Hopf factorization method was initially nsed for solving certain integral equa¬ 
tions [H]. This method leads to Wiener-Hopf equations which can also be recast, by taking 
a Fourier transform, in the form of a Riemann-Hilbert problem [SI [T3] . Although there is a 
well-known uniqueness theory based on the index for Riemann-Hilbert problems, the function 
class usually considered is Holder continuous functions going to one at inhnity, and as will be 
seen, this is not always the right class for applications. We consider the problem of hnding 
positive dehnite solutions to Riemann-Hilbert problems on the real line. This problem has 
a very natural connection with the theory of characteristic functions. This is because, by 
Bochner’s theorem, a positive dehnite function 0 : M"" —)■ C which is continuous at the 
origin and satishes the equation 0(0) = 1, is the characteristic function of some random vari¬ 
able and vice versa. Therefore, positive dehnite functions and characteristic functions are 
nearly synonymous. Furthermore, characteristic functions and probability density functions 
are intimately connected in the sense that probability density functions are inverse Fourier 
transforms of characteristic functions and characteristic functions are Fourier transforms of 
probability density functions. Applications of Wiener-Hopf factorization methods in hnding 
the distribution of certain random variables in statistics dates back to Baxter, Donsker [2] 
and Pecherskii, Rogozin naE]. Since characteristic functions are positive dehnite, we are 
interested in positive dehnite solutions of Riemann-Hilbert problems. This is particularly 
important when the solution cannot be shown to be unique. 

In the usual uniqueness theory for Riemann-Hilbert factorizations, the key condition that 
must be imposed on the function / is that it must go to a nonzero constant at inhnity: this 
condition is needed to insure that the so-called index of the function / is well-dehned and 
stable under small perturbations (see Dehnition [15] and the discussion there). As will be 
seen later, the theory of Levy processes gives natural examples of Riemann-Hilbert problems 
where we seek to factorize a characteristic function as the product of two characteristic 
functions, but the condition of going to a nonzero constant at inhnity is too strong. Thus 
the usual uniqueness theory for Riemann-Hilbert problems becomes problematic. However, 
what is really needed in this application is not so much the uniqueness theory as just the 
existence of a factorization where both factors are positive dehnite, as in the abovementioned 
result. 
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The classical multiplicative homogeneous Riemann-Hilbert problem is to determine a sec- 
tionally analytic function <h with a jump discontinuity on an oriented contour L consist¬ 
ing of a finite number of simple smooth closed or open curves, and with sections $+ and 
<h_ constituting the boundary values on L, which satisfies on L the boundary condition 
G(t) = where the function G, G(t) ^ 0; satisfies the Holder condition every¬ 

where except at a finite number of points where it may have discontinuities of the first 
kind. 

We restrict ourselves to the case where the contour is the real line, because the Bochner 
theory of positive definite functions is only available for a very few contours, the cases of the 
real line and the unit circle are the only ones that appear to be tractable. 

Our main result (Corollary [T^ states that under certain conditions, a function / on the 
real line can be factorized in the form / = $_|_$_ with $_|_ and <h_ are positive definite if 
and only if / is positive definite. 

A signihcant difficulty with the classical singular integral solutions to Wiener-Hopf equa¬ 
tions in particular and Riemann-Hilbert problems in general, is that these solutions often 
involve integrals that are both analytically intractable and difficult to evaluate numerically 
(c.f. [IQ]). In this paper, we develop a method which gives positive definite solutions to 
certain Riemann-Hilbert problems. Furthermore, our method doesn’t suffer from the afore¬ 
mentioned difficulties and gives exact results that can be truncated with a known error term 
to provide accurate numerical approximations directly. 


2. Definitions and Preliminaries 

We dehne positive dehniteness for functions as follows: 

Definition 1. A function on the real line is positive definite if and only if the Fourier 
transform exists everywhere and takes values in the non-negative real numbers. 

We define positive semi-definiteness for matrices as follows: 

Definition 2. A matrix is positive semi-dehnite if and only if it is Hermitian and has 
non-negative eigenvalues. 

It was shown by Bochner [1] that the following Theorem holds: 

Theorem 3. A function f is positive definite if and only if the matrices of the form [f{xi — 
Xj)] are positive semi-definite matrices, for all choices of finite subsets of points {xi} C M. 

Schur has obtained the following interesting property of positive semi-definite matrices 

ng: 

Theorem 4. The elementwise product [ap%] of positive semi-definite matrices [a^] and [%] 
is positive semi-definite. 

Combining these two Theorems, we prove the following: 

Theorem 5. If f{x) is positive definite then exp(/c/) is positive definite (for all scalars 
k > 0). If f{x)^ is positive definite for all A > 0, and f is bounded away from zero, then 
log / differs by a constant from a positive definite function. 

2 


Proof. First of all, we show that if / is positive dehnite, then so is exp{kf). By Theorem [3l 
this is equivalent to showing that if matrices of the form [f{xi—Xj)] are positive semi-dehnite, 
then so are the matrices [exp{kf{xi — Xj))]. But by Theorem 0] we have hrst of all that if 
[f{xi—Xj)] is positive semi-dehnite then the same holds for [/(xj—Xj)™], where m is a positive 
integer. The exponential function has an entire Taylor series expansion exp{kt) = 
where the coefficients are positive, so [exp(fc/(xj —Xj))] can be written as a convergent sum 
of positive semi-dehnite matrices, and is hence positive semi-dehnite. The second claim is 
proven in m p.531], see also [22]. □ 


Finally, we recall the following: 


Theorem 6. Suppose that f has a Fourier transform that is zero outside the interval 
Then, 


/(^) = 


sin f (x — nh) 
f(x-nh) 


hi- 


More generally, if f is continuous, square-integrahle, and is represented as the Fourier trans¬ 
form of a function f in (M) then 


/(^) 


sin f (x — nh) 
^(x — nh) 


is bounded by 



\f{u;)\du. 


The above theorem has a complicated and interesting history, which is reviewed in |2nj . 
There exist a multitude of early versions of the theorem, such as in [2l], or the very early 
paper by de la Vallee Poussin |T3|. The hrst part of the theorem as quoted above comes 
from [19], see also [23], and for the second part see [5] Theorem B]. 

The above Theorem is often used to expand functions whose Fourier transform decays 
at inhnity, rather than actually being zero outside a bounded interval. Thus, for example, 
Rybicki [16] shows by applying the above Theorem that 


F{z) 


lim —= 

h —i-o-i- ^y7^ 


E 

n odd 


exp{ —( 2 ; — 
n 


where F{z) happens to be the complex Dawson’s integral. We now give a Corollary that 
generalizes Rybicki’s observation: 


Corollary 7. Suppose that f G L^(M) is positive definite and continuous. Then, f is 
represented as the Fourier transform of some function in L^(M), and 


fix) 



00 

n=—oo 


sin(7r(x/h — n)) 
'K{x/h — n) 


where the limit is with respect to the supremum norm. 
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Proof. If / is positive definite, then by definition its Fourier transform / is real-valued and 
non-negative everywhere. By the Fourier inversion formula, which holds everywhere since 
/ is continuous, we then have that /(O) = ^ f {cu) du. Since / is non-negative, this 

integral equals the Li-norm of /, and is hnite because the continuity of / insures that /(O) 
is hnite. 

Let us consider the Shannon-Whittaker series expansion for some arbitrary h > 0 : 


OO 

Sh{x) := ^f{nh) 


— OO 


sin f (x — nh) 
^(x — nh) 


We see that if / is continuous, square-integrable, and is represented as the Fourier transform 
of a function / in L^(R) then 

|*^h(x) - /(x)| < - / \f{u)\du. 

Taking the limit as h —)■ 0-1- , we obtain uniform convergence. □ 


We recall the classic form of the Paley-Wiener theorem: 

Theorem 8. Suppose f G L^(M). The following are equivalent: 

(i) The real function f vanishes on R“; 

(a) The Fourier transform f of f extends to a holomorphic function on the upper half¬ 
plane and the L 2 -norms of the functions x i—)■ f^x + iyo) are continuous and uniformly 
hounded for all yo > 0. 


3. RiemANN-Hilbert problem for positive definite functions 


Let L be an oriented contour which consists of a hnite number of simple smooth closed or 
open curves situated in an arbitrary way on the plane. The classical additive (multiplicative) 
homogeneous Riemann-Hilbert problem is to determine a sectionally analytic function $ with 
a jump discontinuity on L and with sections <h+ and <I>_ constituting the boundary values on 
L, which satishes on L the boundary condition <F+(f) — *h-(t) = G{t) ( ‘h+(t) = G{t)^-{t) ) 
where G is a function satisfying the Holder condition everywhere except at a hnite number 
of points where they may have discontinuities of the hrst kind, and G{t) ^ 0 (See [HI section 
44.1]). In the following Theorem, the function / does not need to satisfy the Holder condition. 


Theorem 9. Suppose that f is in L^(R), is positive definite, and continuous. The additive 
homogeneous Riemann-Hilbert problem for f is solved by a decomposition /_ -|- /+ where 
/_i_ has Fourier transform supported in [0, cxo) and /_ has Fourier transform supported in 
(—cx3,0]. The functions in the decomposition are positive definite. We furthermore have 


\ V ^ rf exp(7rz(u;/h-n)) - 1 

t+[w) = hm > t (nh) -;— -r- 

^ h—>-0 ^ ^ 2m(w/h - n) 


and 


f_{w) = ^lim^ 


1 — exp(—27r(t(;/h — n)) 
2ni{w/h — n) 
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Proof. We can decompose uniquely so that / = /_ + /+ where the Fourier transform of /+ 
is the restriction of the Fourier transform of / to R"*", and similarly for /_. If / is positive 
dehnite, then clearly so are /+ and /_. Theorem [8] shows that the holomorphicity conditions 
needed for the additive Riemann-Hilbert problem are satished by /+ and /_. We observe 
that f± is in fact continuous: this is because the hrst part of Corollary [7] shows that the 
Fourier transform / of / is in L^(R). The same is then true for the Fourier transforms of 
/_i_ and /_, and the Riemann-Lebesgue lemma then implies that /+ and /_ are continuous. 
We now consider the boundedness of /+ in the upper half plane. Since /+ is continuous and 
positive dehnite, the hrst part of Corollary [7] shown that f+{x) = du, which 

converges for all x in the upper half plane. But since < 1 for all real uj and all x in the 
upper half plane, we see that |/+(a;)| < |/(i^)| dcj. The series decompositions come from 

a slight generalization of the second part of Corollary [7] to the case of functions supported 
on [0,7r/h] or [—vr/h, 0]. □ 


Proposition 10. Suppose that f : R —> C is positive definite, has Fourier transform 
supported on [—a, a] where 0 < a < tt and 

Im(/(n)) 


lim ita 


= 0 


n —>-±oo n —^±oo Ti 

The function f has a decomposition /_ + /+ where /+ has Fourier transform supported in 
[0, a] and /_ has Fourier transform supported in [—a, 0]. The functions in the decomposition 
are positive definite. We furthermore have 

exp{7ii{w — n)) — 1 


f+{w) = / 


n 


and 


/-H = f{n) 


2'7ii{w — n) 

1 — exp(—— n)) 


2Tii{w — n) 

where the convergence is uniform on compact subsets o/R. 


Proof. The statements on uniform convergence follow from [H Theorem 1]. Since / is positive 
dehnite, f{—x) = f{x) implying that fe{x) = = Re{f{x)) and fo{x) = = 

ilm{f{x)). □ 


Theorem 11. Suppose that f is a function on the real line having the properties: 
i) f e L\R) + C1, 

a) I/I is bounded and bounded away from zero, 

Hi) and f^ is positive definite for all X > 0. 

Then f = where <F+ is positive definite on R and extends to a bounded holomorphic 

function on the upper half plane, and <F_ is positive definite on R and extends to a bounded 
holomorphic function on the lower half plane. 


Proof. By the second and third hypotheses, log(/) dihers by a constant from a positive deh¬ 
nite function (Theorem [5]). It is straightforward to show that the hrst two hypotheses imply 
log(/) is in L^(R) -|-C1. The solution to the additive Riemann-Hilbert problem (Theorem [H]) 
then implies that 

log(/) = g+ + g- + X, 
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where A is a constant, g+ is holomorphic and bonnded in the npper half plane, and g- is 
holoniorphic and bonnded in the lower half plane. Exponentiating, we notice that since the 
exponential fnnction is entire, exp( 5 f_|_ + |) is therefore holomorphic on the npper half plane 
and exp{g_ + is holomorphic on the lower half plane. Bonndedness in these domains 
follows from the bonndedness of g±. Fnrthermore, since g± was positive dehnite on M, the 
hrst part of Theorem [S] implies that := exp( 5 f± + A) is positive dehnite on R. Moreover, 
<h_i_ and <h_ have no zeros in the npper and lower half planes respectively because g^ and g^ 
have no poles in the upper and lower half planes respectively. □ 

Remark 12. The general inhomogeneous problem can then be solved by combining the ho¬ 
mogeneous multiplicative problem with the homogeneous additive problem. 

Corollary 13. If function f satisfies the assumptions of Theorem CD, then f is positive 
definite if and only if there exist positive definite functions and <h_ such that f = <h+<h_. 

Proof. By Theorem flTT]) . if / is positive dehnite then there exist positive dehnite functions 
<h_i_ and <h_ such that / = Its converse holds because product of positive dehnite 

functions is positive dehnite. □ 

In the proof of Theorem (ITT]) , the expression log(/) becomes singular if / vanishes at 
inhnity on the real line. Nevertheless, we can handle this case by perturbing the function 
/ to / = / + e for some e > 0. It is easy to see that the perturbation / also satishes the 
assumptions {ii) and {Hi) of Theorem flTT]) . The following Lemma shows that the process 
of exponentiating and anti-exponentiating which are used in constructing the solution are 
well-behaved with respect to approximation. 

Lemma 14. If x and y are two points in G = {z E C : N < \z\ < M} such that the segment 
joining them is also in G then 

exp(A^)|x — y\ < \ exp(a;) — exp(?/)| < exp(M)|x — y\ 

Proof. By m Theorem 10], exp(a:) — exp(|/) = {x — |/)(Ai exp(r(;i) -|- A 2 exp(tC 2 )) where Ai, 
A 2 > 0, and Ai -|- A 2 = 1 and wi, W 2 are on the segment joining x and y. Therefore, we 
conclude 


exp(iV)|a; — y\ < \ exp(a;) — exp(?/)| < exp(M)|a; — y\ 

□ 

4. An EXAMPLE: DISTRIBUTION OF EXTREMA OF A CLASS OF LEVY PROCESSES 

We thus have a theory that allows us to construct positive dehnite solutions to Riemann- 
Hilbert problems in a very simple and explicit form, well-suited to numerical computation 
when needed. Since positive dehnite functions are, up to a normalization, the same as the 
characteristic functions that appear in statistics, we seek statistical applications. We note 
that in the case of characteristic functions, the Riemann-Lebesgue lemma implies that they 
are in Co(R), and thus the usual uniqueness theory based on the index for Riemann-Hilbert 
problems is not directly applicable. In fact, the index of a function in Co(R) is undehned 
(see [B], however, for some special cases). The usual dehnition of the index is as follows: 
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Definition 15. [Gj Definition 12.1] Let L be a smooth closed contour and G{t) a continuous 
non-vanishing function given on L. By the index of the function G(t) with respect to the 
contour L we understand the increment of its argument, in traversing the curve in the positive 
direction, divided by 27i. 

We will apply the theory discussed in the previous section to the characteristic function of a 
Levy process. We define a Levy process X as in [3l 1. Definition 1]. The function : M —> C 
in 0(A) = E[exp(iAX0] = exp(—^^(A)) where t > 0, is called the characteristic exponent 
of the Levy process X and 0 is the characteristic function which has numerous properties 
such as uniform continuity and positive definiteness. Let r = r(g) be an exponentially 
distributed random time with parameter g > 0 which is independent of the Levy process X. 
The characteristic function of X^. is , see [31 p.l65]. 

Proposition 16. If (p, xf, pi, and a are the characteristic function, characteristic exponent, 
drift, and volatility of a Levy process respectively and a E (0,1), then the following hold: 

i) If t ^ 0 and a ^ 0 then 0 E L^(M) fl L^(M) + Cl and if lim '0(A) = lim A^ then 

A—>-dioo A—>-dioo 

J^gL2(M)nLi(K), 

a) |0| and \ + a\ are bounded and bounded away from zero. 

Hi) (^^)^ H positive definite for all (3 G M'*', 

iv) lndex(0 + « + !) = Index(^^ + a) = lndex(0 + a) = 0. 

Proof, i) By the Levy-Khintchine formula, 

/ OO 1*1 POO 

(1 — cos{xX))i'{dx) H—cr^A^ + *( / xXiy{dx) — / sin(a;A)z/((ia;) — /lA) 

OO 2 J _l J — oo 

where z/ is a Levy measure. Since f(cos(a;A) — 1) < 0 we conclude Re(—f0(A)) < — 
Therefore, 



|0(A)|dA < 



-ta^X^)dX 



^ta‘^X‘^)dX 


hm - — 

x^Qo 2a\/t 



lim 

A— y —OO 


2a\/t 




Similarly, 


|0(A)|^(iA= / exp(—2f Re(0(A)))(iA < 


TT 


a 


Vt 


It is also well-known that a uniformly continuous absolutely integrable function vanishes 
at infinity. If [3l Proposition 2 (i)] holds then 0 behaves like A^ at infinity implying that 
goes to zero at infinity like and consequently belongs to L^(M) fl L^(M). 
ii) Since 0 < |0| < 1, 0 is bounded and bounded away from zero. Since Re(0) > 0, we 
conclude from [3l Proposition 2] that -|- a is bounded and bounded away from zero. 

Hi) Since is a characteristic function, it is positive definite by definition. Let 0 > 0. 

From the definition of the gamma function we have 
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\P - 




'q + ^{X)' 


r(/9) 


exp(- 


,g +V’(A), 


x)dx = 




exp(-a;) 

- exp{-{-)'ilj{X))dx 


The function exp(—(|)'^(A)) is positive definite with respect to A and is positive. 

Therefore, Y is positive dehnite for all /9 > 0. 

iv) It is well known that the modulus of the characteristic function is less than or equal to one 
and since the modulus is also equal to exp(—f Re('0)) and f > 0, we conclude that Re('0) > 0. 
This fact together with [3l Proposition 2] indicates that the functions 0 + a + l,'0 + Q!, and 
+ a are all non-vanishing. Therefore, Re('0 + a) > a implying that index oi 'ip + a is 
zero. Similarly, the curve traced out by 0 -|- a -|- 1 = exp(—f Re('^)) exp(—R Im('^)) + a + 1 
resides in the right half plane and doesn’t enclose the origin. We furthermore observe that 


Re( 


g +'?/’ 


+ “) — ( 


|g +V’l 


f + a + 


q + 'ip\- 


Re(V’) 


Since Re(^:^ + «) > (Y + a, we conclude that Index(^^ -|- a) = 0 


■ k+bl 


q+Tp 


□ 


Part iv of the above proposition is just included for its intrinsic interest: we do not 
actually use the index of T -|- a for anything. We next solve the multiplicative Riemann- 
Hilbert problem for the function of the above proposition. This function goes to zero 
at inhnity (quadratically) and therefore does not have an index. Thus, we cannot argue 
that there exists a unique and positive dehnite solution(s) to the Riemann-Hilbert problem, 
but we can use our techniques to at least hnd a positive dehnite solution. The supremum 
and inhmum of the Levy process where r = T{q) is exponentially distributed are dehned 
to be Mg = sup{Xs : s < r(g)} and Ig = mfjXs : s < r(g)} respectively. It is shown in 
[3[ VI.2.Theorem 5] how to solve for the supremum and inhmum in terms of a Riemann- 
Hilbert factorization problem, and evidently the factors must be positive dehnite to be valid 
solutions. The next theorem thus gives the distributions of the extrema. This theorem is 
more general and it is not limited to particular classes of Levy processes such as a-stable 
processes considered in [7]. 

Theorem 17. The distributions fMq and fi^ of the extrema Mg and Ig of the Levy process 
Xj. where r = r(g) is exponentially distributed are given as follows: 

i) fMq{X) = 'ifq+iX) where iiq+{X) = exp{{iX 1)5'+(A) - f), and 


g+{X) = lim g{nh) 


exp7ri(A/h — n) — 1 
h — 5.0 ^^ ' ‘IniiXjh — ri) 

n=—oo ^ ' 

a) fiq{X) = 'ipq_{X) where 'ipq_{X) = exp{{iX + l)g_(A) — |), and 

1 — exp —i7i{X/h — n) 


g_{X) = ^lim^ Y dinh) 


2m{X/h — n) 


The function g is defined by g{X) = ^+in(<?)-M 9 +h(A)) ^ 















Proof. By Proposition flT6|) {in), is positive definite for all /? > 0 implying that 

In(q') — ln(q' + ipiX)) differs by a constant from a positive definite function. Therefore, there 
exists a constant c such that c + ln(g) — ln(g + 'if{X)) is positive dehnite. Although this 
function is not in L^(R), we will multiply it by another positive dehnite function so that 
the product function belongs to A suitable function is which is positive dehnite 

because its inverse Fourier transform is non-negative. According to Theorem ([9]), there exist 
positive dehnite functions g+ and g- such that g{X) = = g+{X) + g-{X). By 

exponentiating we have: 

= exp((iA + 1)^+(A) - |) exp((iA + l)^-(A) - |) 

Since multiplication by {iX + 1) corresponds to applying a diherential operator to the 
Fourier transform, and since diherential operators do not increase support, it follows that 
{iX + l) 5 f+(A) and {iX + l)g-{X) are positive dehnite. Thus we obtain a factorization 
'f’q+{X)'ipq_{X) of If is shown in [3l VI.2.Theorem 5] that the hrst factor is the char¬ 

acteristic function of Mq and the second factor is the characteristic function of Iq. The 
distributions are given by Fourier transform of the characteristic functions. 

□ 

Theorem 18. Let r = T{q) he a geometrically distributed random time with parameter 
0 < q < 1 and let 0(A) = exp(—-(/^(A)) be the characteristic function of the Levy process Xi 
such that lim '0(A) = lim A^. Then, the characteristic function of Xr which is of the 

A—>-ioo A—>-ioo 

form /(A) = i^q^(^x) following properties: 

i) f e L\R) + C1, 

a) I/I is hounded and bounded away from zero, 

Hi) f^ is positive definite for all (3 > 0. 

Proof, i) The function / decays to 1 — g like an exponential. Therefore, / G L^(M) -|- Cl. 
a) Since 0 < |0(A)| < 1, it follows that 1 — g < |/(A)| < 1 for all A G M. 

Hi) By the binomial series, we have 


1 

(1-#(A))/^ 




k=0 



k 


Since ^ j = is positive, g^ and (1 — g)^ are positive for all 

(3 > 0 and (0(A))^ = exp(—fc0(A)) is positive dehnite, we conclude that is positive dehnite 
for all 0 > 0. □ 


The following corollary follows from the above Theorem together with Theorem flTT]) : 

Corollary 19. Let r = r(g) be a geometrically distributed random time with parameter 
0 < g < 1 and let 0(A) = exp(—0(A)) be the characteristic function of the Levy process Xi 
such that ^\im 0(A) = ^lim A^. Then, the characteristic function ofX^, i.e. /(A) = 

has a factorization of the form /(A) = where <F+ is positive definite on M and extends 

to a bounded holomorphic function on the upper half plane, and <h_ is positive definite on M 
and extends to a bounded holomorphic function on the lower half plane. 
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Proposition 20. Let {L~^,H) be the ladder process of a given Levy process X with charac¬ 
teristic function (f){a,(3) = E[exp(iaL“^(t) -\-i(3H(t))] = exp{—K{a, (3)t) where t > 0 and the 
bivariate Laplace exponent K{a,/3) is of the following form [3l VI.2.Corollary 10].- 

dt f (e-* - e dx) 

J [0,cxd) 

Then, the following hold: 

i) \(j)\ is bounded and bounded away from zero, 
a) (fP is positive definite for all ■j > 0. 

Proof, i) Since 0 < |0| < 1, 0 is bounded and bounded away from zero. 

ii) Since 0 is a characteristic function, it is positive definite. Hence, K,{a,(3) is negative 
definite. Since 7 > 0 , 'yK{a,/3) is negative definite which implies that (p is positive dehnite. 

□ 


K{a, (3) = kexp ^ J 
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